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The non-modal self-heating mechanism driven by the velocity shear in kinematically complex 
magnetohydrodynamic (MHD) plasma flows is considered. The study is based on the full set of 
MHD equations including dissipative terms. The equations are linearized and unstable modes in 
the flow are looked for. Two different cases are specified and studied: (a) the instability related to 
an exponential evolution of the wave vector; and (b) the parametric instability, which takes place 
when the components of the wave vector evolve in time periodically. By examining the dissipative 
terms, it is shown that the self-heating rate provided by viscous damping is of the same order of 
magnitude as that due to the magnetic resistivity. It is found that the heating efficiency of the 
exponential instability is higher than that of the parametric instability. 

PACS numbers: 98.62. Nx, 98.54.Cm, 94.20.wf, 95.30.Qd 



I. INTRODUCTION 

It is widely known that astrophysical, terrestrial and 
laboratory plasma flows are most often characterized by 
spatially inhomogeneous velocities, constituting so-called 
shear flows. In many cases, these flows are kinematically 
complex, multidimensional and sometimes they are also 
relativistic. Among the different kinds of plasma shear 
flows many exhibit motion profiles with highly nontrivial 
velocity fields. One prominent class of such flows con- 
tains the swirling astrophysical flows that until now have 
received only little theoretical attention. In particular, 
many galaxies reveal relativistic jets [lj, some of which 
are characterized by helical motion patterns @, Q. A 
similar situation is met in pulsars, where morphological 
and spectroscopic studies make strong jet-like structures 
evident 0,@. Even closer to us, namely in the solar at- 
mosphere, the existence of giant macro-spicules featuring 
swirling, tornado-like plasma motions has been reported 

!• 

The latter class of jet-like large-scale velocity patterns 
in the atmosphere of the Sun deserves some special at- 
tention. In 1991, the Yohkoh Solar Observatory detected 
coronal jets (by means of the Soft X-Ray Telescope Q) 
which were assumed to be related to evaporation flows 
Q. More recently, it has been reported on the 3 He- 
rich solar energetic particle event [9[ observed on 2006 
November 18 th by the X-Ray Telescope on board the 
Hinode satellite There is also some indirect obser- 
vational evidence for the presence of more complicated 
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velocity patterns. In particular, blue shifted and red 
shifted emission on either side of the axis of a macro- 
spicule has been reported @. These observations were 
interpreted as indicating at the presence of a rotation 
within these tall pillars of hot and magnetized, swirling 
plasma flows. Another interesting class of astronomical 
objects with an observational evidence of swirling flows is 
the Herbig-Haro jets. It was shown that the young stellar 
object HH 212 exhibits strong indications for the pres- 
ence of rotation in the jet [III . A similar helical structure 
is observed in DG Tau fl2l | . 

Therefore, the study of kinematically complex shear 
flows with a non-trivial geometry is interesting and neces- 
sary because flow patterns with nontrivial velocity fields 
appear to be quite common in astrophysics. Besides, it 
seems reasonable to expect that the kinematic complex- 
ity could considerably influence the physical processes 
within these flows and lead to various detectable obser- 
vational appearances in the related astronomical objects. 

Recently, it was fully realized that collective phenom- 
ena in shear flows are characterized by so-called non- 
modal processes, related with non- normality These 
processes are fairl y w ell-understood for a variety of flows, 
see, for instance, |T3 - [r7| . In particular, in [li| [I?} shear 
flow instabilities are studied in helical flows. The problem 
was considered both for the incompressible [1 61 ) and com- 
pressible [ljj cases and it has been shown that shear flows 
efficiently amplify the fast/slow magnetosonic waves (in 
the case of compressible flows) and the Alfven waves (in 
both cases). It was also argued that these instabilities 
may lead to so-called 'non-modal self-heating'. This pro- 
cess involves three major stages: (a) First, the waves are 
excited spontaneously within the shear flow; (b) second, 
the excited modes undergo a non-modal amplification, 
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extracting a part of the equilibrium flow's kinetic energy 
and (c) third, in the final stage, the non-modally am- 
plified waves get dissipated via viscous damping and/or 
magnetic resistive diffusion and give their energy back to 
the background flow in the form of heat. 

This mechanism was originally described in [l8| for 
the case of a hydrodynamic flow and it has been ar- 
gued that self-heating may play an important role also in 
the dynamics of kinematically complex magnetized flows. 
In [l9j it has been shown that efficient non-modal self- 
heating can indeed take place in magnetized media as 
well. A direct application to the solar atmosphere was 
proposed in (2(| where it was concluded that self-heating 
could serve as an efficient mechanism for providing, or 
at least contributing to, the heat production in the lower 
corona. Recently, we have studied the efficiency of non- 
modal self-heating by acoustic wave perturbations in the 
solar chromosphere [21] . It has been shown that ampli- 
fied waves get damped due to the presence of viscous 
dissipation and indeed cause an energy transfer back to 
the background flow in the form of heat. 

In most of the above cited papers, the problem of self- 
heating was considered only for shear flows with rela- 
tively simple flow geometry, while the flows occurring in 
astrophysical plasmas, as well as in many terrestrial and 
laboratory situations, can be significantly more complex. 
Therefore, it is quite timely and important to examine 
kinematically more complicated cases and verify if the 
mechanism is also efficient when such complex flow pat- 
terns occur. For this purpose, in the present paper we 
study the possibility of non-modal self-heating in fully 
compressible, kinematically complex MHD flows taking 
into account both viscosity and magnetic resistivity. 

The paper is arranged in the following way. In sec- 
tion II, we formulate the model for the study of non- 
modal shear flow instabilities for MHD flows including 
resistive terms and we derive the equations governing the 
self-heating process. In section III, we apply our model to 
the fully compressible flows examining the high plasma- 
f3 case and we then investigate numerically the efficiency 
of the self-heating mechanism. In section IV, we discuss 
the obtained results and outline the directions of further, 
more detailed, quantitative and applicative studies. 



II. MAIN CONSIDERATION 

In order to study the self-heating mechanism in mag- 
netized plasma flows, we start with the standard set of 
polytropic MHD equations consisting of the equation of 
mass conservation: 

D t p + pV ■ V = 0, (1) 
the momentum conservation equation: 
1 B 

D t V = — VP x (V x B) + vAV, (2) 

p Anp 



the induction equation: 

D t B = (B • V)V — B(V • V) + r?AB, (3) 

Gauss's law for magnetism: 

V • B = 0, (4) 

and the polytropic equation of state: 

P = Cp 1+1/n , (5) 

where Dt = dt + (V • V) denotes the convective deriva- 
tive, p is the density, P stands for the pressure, V is the 
velocity field, B denotes the magnetic field, v stands for 
the coefficient of kinematic viscosity, r\ is the magnetic 
resistivity coefficient, and C = const and n are the poly- 
tropic constant and the polytropic index, respectively. 
Note that in the framework of the present paper we con- 
sider a common stellar model, viz. the so-called "Poly- 
tropic Star Model" in which the pressure depends only 
on the density according to the polytropic equation of 
state, Eq. (5) 

The equilibrium state of the model is specified by a 
homogeneous MHD plasma (po = const), embedded in a 
uniform, vertical magnetic field (Bo = [0, 0, Bq = const]). 
For analyzing the full set of MHD equations, we first 
linearize them, perturbing all physical quantities about 
the equilibrium state: 

p = p a +p, (6) 
V = U + v, (7) 
B = B + B', (8) 

and 

P = P + P'. (9) 

It is assumed that the perturbations p' , v, B' and P' are 
small compared to the corresponding equilibrium quan- 
tities, viz. po, U, Bo and Pq, respectively. By substitut- 
ing Eqs. (OH]) into Eqs. ([T][5]), and then linearizing (i.e. 
keeping only the first-order terms in the perturbed quan- 
tities), one can derive the following set of equations for 
the perturbed quantities: 

V t D + V • v = 0, (10) 
X> t v+(vV)U = -C 2 s VD + C A [d z h-Vb z ] + v&v, (11) 
V t h = (b ■ V)U + d z w + e z (V • v) + ?/Ab, (12) 

V-b = 0, (13) 

w here V t = d t + (U • V), D = p'/p , b = B'/B , C s = 
y/ dPo I dp a denotes the sound speed and Ca = Bq/ y/4irp a 
is the Alfven speed. 
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Following the method developed in [23j, we assume 
that U is a spatially inhomogeneous vector velocity field. 
We expand the velocity in Taylor series around a point 
A(xo,yo, zq), preserving only the linear terms: 



\J = V(A) + ^ 



dU{A) 



(14) 



where i = 1,2,3 and Xi = [x, y, z). 

Mahajan and Rogava have shown that by applying the 
following ansatz 



F(x,y,z,t) = F(i)e 01 ^ 2 



Y,Ki(t)xi 



3 

J2 u i( A ) J Ki(t)dt, 



(15) 



(16) 



(17) 



to Eqs. (|10lllip , these equations transform to a set of ordi- 
nary differential equations [23j j . This reduces the problem 
mathematically to the study of an initial value problem. 
By Ui(A) we denote the equilibrium velocity components 
and the Ki(t) represent the wave number vector compo- 
nents satisfying the following set of differential equations 



<9 t K + S • K = 0, 



(18) 



where S T is the transposed shear matrix. This shear 
matrix S is defined by: 



Uy^x &y,y Uy,z 



(19) 



with Ui t k = dUi/dxk- 

Working with the background velocity field U(r) = 
[0, rfl(r), U(r)], the shear matrix obtains the following 
form [la, Il7ll 



£ At 
A 2 -S 
Ci C 2 



(20) 



By taking Eqs. 



into account, the set of 



Eqs. (|10H13p reduces to the following dimensionless form 



dW=k-u, 



(21) 



u (1) + s • u = -ke 2 d + h - kh z - 77fc 2 u, (22) 



h (1) =s-h-u-e z (k-u)-^fc 2 h, (23) 



k-h = 0, 



• k = 0, 



with 



(24) 



(25) 



(26) 



where d = iD, u = v/C A , h = ih, e = C s /Ca, 
a x , 2 = A lt2 /K Z C A , a = Y,/K Z C A , n, 2 = C h2 /K Z C A , 
k = K/K z , V = vK z jC A and 77 = r,K z /C A . By we 
denote dF/dr and r = K z C A t is the dimensionless time. 

For estimating the efficiency of the self-heating mech- 
anism, it is useful to define the total energy of the per- 
turbations: 




E,, 



-&kin ~~r~ -&r> 



E r 



u 

T 



h_2 

2 



e 2 d 2 



(27) 



which consists of the kinetic, magnetic and compressional 
energies, respectively. The evolutionary equation for the 
total energy has the following form: 



E 



(i) _ 



-u(s • u) + h(s • h) - 2h z (k ■ u)- 



—1 2 2 
vk u 



rjk 2 h 2 . (28) 



As it is clear from Eq. (f28|). the last two (viscosity and 
magnetic resistivity) terms are responsible for the energy 
dissipation and thus for the final stage of the self-heating 
process. By means of these two terms one can define the 
so-called 'self- heating rate', originally introduced in [l8[ 
only for purely hydrodynamic flows: 

* (T) = E^fij [ [^VJuV) +^kV)h 2 (r')] dr'. 

(29) 

III. DISCUSSION 

In this section, we study the efficiency of the self- 
heating mechanism for different values of parameters, 
taking into account both the viscous damping and mag- 
netic resistivity terms. Generally speaking, for the non- 
modal self- heating process to be efficient, the dissipative 
terms must not be too small in order to provide substan- 
tial heating. However, at the same time the damping 
should not be "too strong" either: non-modally ampli- 
fied modes have to have enough time for growing and 
extracting the equilibrium flow energy before damping 
occurs. As it is clear from Eqs. (J5J3]), the effects of dissi- 
pation become important only for relatively small length 
scales. 

Mathematically the non-modal instability is described 
by the system of equations (|2HI25[) . where Eq. (|2"51 plays 
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FIG. 1. The temporal behavior of h x (r) and the energy 
E n (r) = E to t(r) / E to t{Q) normalized by its initial value for the 
following set of parameters: a — 0, 01 = — an = — 1, ri = 0.01, 
T2 = 0.9406, e = 0.1, k x o = fc^o = 10, u x o = u y o = u z o = 0, 
h x o = 0.01, h y o = 0, do = 0, V = rj — 0. Note that the range 
of r2, where the evolution of the modes is parametrically un- 
stable, is very narrow: [0.94; 0.941]. 



an important role. In particular, in plane-parallel flows 
k exhibits a linear time evolution [24j, whereas for mul- 
tidimensional cases (such as for helical flows), the com- 
ponents of the wave vector may evolve in time either 
exponentially or periodically [lot Il7| . Linearly and ex- 
ponentially evolving wave vectors lead to so-called 'usual' 
instabilities (with T 2 = a 2 + a\a2 > 0). Parametri- 
cally unstable modes deserve a special interest. These 
modes appear for periodically evolving k's (i.e. T 2 < 0) 
for certain values of the parameters. In Fig. [T] we show 
the temporal behavior of h x (r) and the energy E n {r) = 
Etot(j)l Etot(0) normalized by its initial value. The con- 
sidered set of parameters is: a = 0, a% — —02 = — 1, 
ri = 0.01, r 2 = 0.9406, e = 0.1, k xQ = k y0 = 10, 
Uxo = u y o = u zQ = 0, h x0 = 0.01, h y0 = 0, d Q = 0, 
77 = rj = 0. As it is clear from the plots, despite T 2 < 0, 
the system exhibits unstable modes. The instability takes 
place only for very narrow ranges of parameters. In 
particular, one can straightforwardly check that the in- 
stability disappears for values of r 2 outside the interval 
[0.94; 0.941]. 

The dissipative factors will inevitably lead to the self- 
heating effect. In particular, Fig. [5] shows the temporal 
evolution of the self-heating rate for two different cases: 
on the top panel we examine V = 1.25 x 10 -4 , rj = 
and on the bottom panel the result is shown for v = 
and rj = 1.25 x 10~ 4 . All the other parameters are the 
same as in Fig. [TJ As we see from the plots, there is 
no principal difference between viscosity-dominated and 
magnetic-resistivity-dominated cases, i.e. the saturated 
values of the self-heating rates are of the same order of 
magnitudes. In both cases, the energy converted into 
heat exceeds the initial perturbation energy about 10 
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FIG. 2. The temporal evolution of the self-heating rate for 
two different cases. Top panel: V = 1.25 x 10~ 4 , rj = 0. 
Bottom panel: V — 0, rj — 1.25 x 10~ 4 . The set of parameters 
is a = 0, oi = -a 2 = -1, n = 0.01, r 2 = 0.9406, e = 0.1, 
k x o = k v o = 10, Uxo = u y o = Uzo = 0, hxo = 0.01, h y0 = 0, 
d = 0. 

times. This means that the excited non-modal waves 
extract energy from the mean flow and this energy, via 
the agency of the dissipative terms, is converted into heat 
and given back to the flow. As we see from the results, 
the efficiency of the self-heating is quite high. 

It is worth noting that for the considered mechanism 
to be efficient, several conditions must be satisfied. At 
the one hand, the instability must be robust enough to 
provide an amplification of the waves and, additionally, 
the dissipation process must not be too efficient in order 
to avoid damping of the excited modes before they grow 
substantially enough. On the other hand, if the damping 
is too weak the amplified waves fail to give back their 
energy back to the flow. Therefore, one may expect that 
self-heating will be efficient for a limited range of mod- 
eratly dissipative systems. 

On the top panel of Fig. El we show the temporal evolu- 
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FIG. 3. The temporal evolution of the self-heating rate for 
three different cases. Top panel: 77 = 1.2 x 10" 5 (solid line) 
and V = 1.24 x 10~ 5 (dashed line). Bottom panel: V = 1.22 x 
10 -J . The set of parameters is a = 0, ai — —ai = — 1, 
n = 0.01, r 2 = 0.9406, e = 0.1, k x0 = k y0 = 10, it^o = %o = 
u-zo = 0, /i^o = 0.01, h y o = 0, do = 0. 



tion of the self-heating rate for 77 — 1.2 x 10~ 5 (solid line) 
and 17 = 1.24 x 10 -5 (dashed line). The other parameters 
are again the same as in Fig. [TJ As it is evident from the 
plots, for the case presented by the solid line, the viscos- 
ity is not strong enough and non-modaly amplified modes 
avoid any significant viscous dissipation. In general, the 
temporal evolution of the self-heating rate is very sensi- 
tive to viscosity. By slightly increasing the value of the 
viscosity, the situation might change drastically. By the 
dashed line we show the result for V = 1.24 x 10 -5 , and 
it is clear that unlike the previous case, more dominant 
viscous terms here ensure stronger self- heating. The ef- 
ficiency of the instability seems significantly decreased 
and \&(t) tends to become saturated (see Fig. [2]). This 
tendency implies that for a certain value of the viscosity 
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FIG. 4. The temporal behavior of u x (t) and E n (r). The 
considered set of parameters is a = 0, ai = a2 = 1, r\ = 
r% = 0, e = 0.1, k x0 = 1, kyo = 0, u x0 = 0, u y0 = u z0 = 0, 
hxo = h y o = 0, d = 0, 77 = 77 = 0. 



the effects of non-modal instability and dissipation will 
balance each other. This particular case is presented on 
the bottom panel of Fig. [HJ where the value of the vis- 
cosity equals 1.22 x 10 -5 and, as we see from the figure, 
the self-heating rate shows a linear dependence on time. 

We have already noticed that apart from the paramet- 
ric instability, the system under consideration exhibits 
also 'usual' instabilities (r 2 > 0), characterized by ex- 
ponentially evolving wave number vector. Normally, in 
these cases the waves are more unstable than the para- 
metrically unstable modes. In particular, in Fig. 2] we 
show the temporal behavior of u x (t) and E n {r). The 
considered set of parameters is: a = 0, a\ = 02 = 1, 
n = r 2 = 0, e = 0.1, k x0 = 1, k y0 = 0, u x0 = 0, 
u y o = u z0 = 0, h x o = h y o = 0, d — 0, 77 = rj = 0. 
As we see from the figure, the energy perturbation and 
u x significantly increase until r = 40. These results can 
be explained as follows: since T 2 > 0, the wave vector 
components evolve with time exponentially (unlike the 
cases presented by figures [T][3]) and correspondingly, the 
system undergoes a strong instability, the growth rate 
of which exceeds that of the parametric instability. For 
the same reason, the efficiency of the self-heating mecha- 
nism for the unstable non-modal modes shown in Fig. [4] 
must be higher with respect to the self-heating rate for 
parametrically unstable waves. 

In Fig. [SJ we show the temporal evolution of ^(r) and 
E n (r). The considered set of parameters is the same 
as in the previous figure, except for 77 = 1 x 10~ 5 . It 
is clear from the plots, that initially the instability rate 
exceeds that of the self-heating (negligible role of the dis- 
sipative factors), and therefore, the perturbation energy 
increases. In the course of time, due to the exponen- 
tial growth of the wave vectors' components, the effective 
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FIG. 5. The temporal behavior of *(t) and E n (r). The 
considered set of parameters is a = 0, ai = ai — 1, n = 
r 2 = 0, e = 0.1, k xQ = 1, k y0 = 0, u x0 = 0, M a0 = u z0 = 0, 
fe»o = h y0 = 0, do = 0, V = 1 x 10~ 5 , rj = 0. 



length scale, A ~ 1/fc decreases, leading to an increasing 
role of the viscous terms. As a result, the waves be- 
come getting damped, the perturbation energy reaches 
its maximum value, decreases and finally vanishes (see 
Fig. [5]). This in turn, leads to a self- heating rate with a 
saturated value, ~ 270 — 280, which is much higher than 
^oo for the previous cases. 

We study compressible MHD flows and all results are 
obtained for e = 0.1. We did not discuss other regimes 
(like e ~ l,e 3> 1), because as it turned out, there is 
no qualitative difference between the results obtained for 
this and those regimes. 



IV. CONCLUSIONS 



dissipative terms, it has been shown that for proper val- 
ues of viscosity and magnetic resistivity the self-heating 
rate may be of the order of 10 (see Fig. [5]). The self- 
heating mechanism has also been studied former class of 
instabilities and it was found out that in this case the 
efficiency of self-heating is considerably higher than for 
the case of parametrically unstable modes (see Fig. [SJ . 
This is directly related to the fact that in the time inter- 
val from the initial appearance of these modes until the 
moment, when the dissipative terms become important, 
exponentially evolving modes manage to amplify more 
efficiently than periodically evolving (parametric) modes 
(see Fig. [4]). Consequently, the energy given back to the 
background flow in the form of heat is larger in the for- 
mer case and the resulting self-heating is more efficient 
too. 

The present results show high efficiency of non-modal 
sclf-hcating in kinematically complex MHD shear flows, 
which may be of astrophysical relevance for different 
classes of astronomical objects. In particular, the so- 
lar atmosphere hosts various kinds of jet-like structures 
possibly characterized by helical motion and swirling ve- 
locity patterns On the other hand, it is worth noting 
that in the solar atmosphere there is the long-standing 
problem of chromospheric and coronal heating [2^, [2(| . 
A similar problem of unclear origin of heating source ex- 
ists in the young stellar objects, particularly in Herbig- 
Haro jets. Therefore, it is of great importance to check 
whether the non-modal self-heating could be one of the 
mechanisms of plasma heating in these astrophysical situ- 
ations. Obviously it remains to see on a quantitative level 
if the self-heating in MHD flows with kinematic complex- 
ity may account for the actually observed heating in con- 
crete, realistic cases of astrophysical interest. Similarly, 
the results obtained in this paper can be relevamt and 
important for terrestrial (laboratory) plasma MHD dis- 
sipative and resistive flows where the presence of velocity 
shear and kinematic complexity is known. 



We have considered velocity shear induced instabili- 
ties in viscous and resistive, kinematically complex MHD 
flows. By linearizing the system of equations, two differ- 
ent kinds of instabilities have been disclosed and stud- 
ied: (a) the instability related to exponentially evolving 
wave number vectors; and (b) the parametric instabil- 
ity related to wave number vectors with periodic time 
dependence. The latter instability works only in a rela- 
tively narrow area of parametric space. In particular, it 
has been found that in 3D shear flows one can find spe- 
cific ranges of parameters, where the system undergoes 
a very unstable regime (see Fig. [T]). By considering the 
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